Roton-type mode softening in a quantum gas with cavity-mediated long-range 

interactions 



(N 

o 

oo 



i 

c 

c3 



i 

C 

o 
o 



> 

(N 
(N 

m 

en 

o 

> 

X 




supersolid 
phase 




R. Mottl, F. Brennecke, K. Baumann, R. Landig, T. Donner, and T. EsslingeiQ 
Institute for Quantum Electronics, ETH Zurich, 8093 Zurich, Switzerland 

Long-range interactions in quantum gases are predicted to give rise to an excitation spectrum 
of roton character, similar to that observed in superfluid helium. We investigate the excitation 
spectrum of a Bose- Einstein condensate with cavity-mediated long-range interactions, which couple 
all particles to each other. Increasing the strength of the interaction leads to a softening of an 
excitation mode at a finite momentum, preceding a superfluid to supersolid phase transition. We 
study the mode softening spectroscopically across the phase transition using a variant of Bragg 
spectroscopy. The measured spectrum is in very good agreement with ab initio calculations and, 
at the phase transition, a diverging susceptibility is observed. The work paves the way towards 
quantum simulation of long-range interacting many-body systems. 
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Experiments with ultracold gases have succeeded in 
realizing a variety of quantum many-body phases by ex- 
ploiting the tunability of short-range atom-atom inter- 
actions pp. Creating quantum gases with long-range in- 
teractions 2 - 5 is motivated by the prospect of observing 
previously unexplored phenomena and phases [H] . In par- 
ticular, long-range interactions in a Bose-Einstein con- 
densate (BEC) have been predicted [7HTU] to give rise 
to an excitation spectrum similar to the roton spectrum 
observed in superfluid helium |11) . A roton spectrum 
emerges from density-density correlations which can be 
induced by short-range van-der-Waals interactions as in 
liquid helium or by momentum-dependent long-range in- 
teractions in dilute quantum gases [T . Such a mode soft- 
ening at finite momentum has been proposed as a possible 
route to a supersolid phase [12"] . 

Important steps towards realizing quantum gases with 
long-range interactions have been achieved by cooling 
atomic species with large magnetic dipole moment to 
quantum degeneracy [5J [5] and by the efforts to prepare 
ultracold polar molecules in their ground state U ■ Ef- 
fects of long-range interactions in quantum gases have 
been observed in the anisotropic expansion, in collective 
as well as Bloch oscillations, and in the stability of inter- 
acting dipolar BEC [TSHTo] . However, a roton-type mode 
softening requires very strong long-range interactions and 
is challenging to observe in quantum gases. 

A different approach to long-range interactions in cold 
gases makes use of the radiative coupling between elec- 
tric dipoles induced by off-resonant laser light, which has 
been shown theoretically to cause a roton-like minimum 
in the energy spectrum of a BEC [HI HO] ■ A large enhance- 
ment of such radiative coupling is achieved by placing the 
dipoles into an optical resonator, which leads to strong 
global atom-atom interactions [17l [HI [2Ql [33] . 

We create long-range interactions in a BEC by plac- 
ing it into an optical high-finesse cavity and irradiating 
the atoms with a transverse pump laser, which is far de- 
tuned from the atomic transition frequency [3TJ [351 131)] 
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FIG. 1. Experimental scheme and mode softening. (A) 
A 87 Rb BEC inside a Fabry-Perot resonator is transversally 
illuminated by a far red-detuned standing- wave laser field. 
In a quantized picture, atoms off-resonantly scatter pho- 
tons from the pump field into a close-detuned TEMoo cav- 
ity mode and back at rate NV/h (see text), creating and 
annihilating pairs of atoms in the superposition of momenta 
(Px,Pz) = (±ftfc, ±hk) (see close-up displaying one of four pos- 
sible processes). This results in global interactions between 
all N atoms. The interaction strength V is controlled via the 
power of the transverse laser field. (B) The cavity-mediated 
atom-atom interaction causes a softening of a collective exci- 
tation mode at momenta (±hk, izhk), and a diverging suscep- 
tibility (blue shade) at a critical interaction strength (dashed 
line). 



(Fig. 1 A) . The induced electric dipoles of the atoms oscil- 
late at the pump laser frequency and collectively couple 
to a single cavity mode, which is detuned from the pump 
frequency by a few cavity linewidths. In turn, the in- 
duced cavity field leads to an ac Stark shift in the atoms. 
The resulting atom-atom interaction extends over the en- 
tire atomic cloud and is tunable in strength. For a critical 
strength of the long-range interactions, the BEC under- 
goes a quantum phase transition to a supersolid phase 
with checkerboard density order [3S] , as observed in [32] . 
At the phase transition a discrete spatial symmetry pro- 
vided by the cavity mode structure is broken, thus es- 
tablishing nontrivial diagonal long-range order. We have 



developed a method to measure the excitation spectrum 
across this phase transition. The method combines Bragg 
spectroscopy [24] and cavity-enhanced Bragg scattering, 
and we use it to observe a notable change in the disper- 
sion relation while crossing the phase transition. This 
is complemented by a diverging susceptibility, in agree- 
ment with the second-order nature of the phase transition 
(Fig. IB). 

The cavity-mediated long-range interaction is de- 
scribed by the Hamiltonian 

ff aa = J#t( r )^t( r ')v( r ,r')*(r')*(r)d 3 rd 3 r' with 
atomic field operator ^(r). This Hamiltonian is obtained 
from the dispersive atom-light interaction Hamiltonian 
by adiabatically eliminating the fast cavity field dynam- 
ics 25, 33J . The interaction potential has the form 

V(r, r') = V cos(kx) cos(kz) cos(kx') cos(kz') (1) 

and describes the ac Stark shift experienced by an atom 
at position r as a result of the cavity field induced by 
a second atom at position r'. As a consequence of the 
interference between the mediating cavity field and the 
transverse pump field, the spatial dependence of the in- 
teraction is determined by the corresponding mode func- 
tions cos(kx) and cos(fcz), where k = 2ir/\ denotes the 
optical wavevector (Fig. 1A). The interaction strength 
V = hrj 2 /A c depends on the two-photon Rabi frequency 
r\ which can be tuned experimentally via the transverse 
pump power P |25j . The sign of V is determined by 
the detuning A c = lo p — ui c between the transverse pump 
laser frequency lu p and the dispersively shifted cavity res- 
onance u) c . In the experiment |A C | is large compared to 
the cavity linewidth 2k. 

For negative V, the cavity-mediated interaction in- 
duces density correlations in the atomic cloud with spa- 
tial periodicity of A along the pump and cavity direction. 
In momentum space, this corresponds to the creation and 
annihilation of pairs of correlated atoms in the momen- 
tum mode |e), which is the symmetric superposition of 
the four momentum states \p x ,Pz) — I ± hk,±hk). For 
a macroscopically populated zero-momentum mode con- 
taining N atoms, this is described by the Hamiltonian 
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FIG. 2. Probing the excitation spectrum in the normal phase 
(P = 0.65P cr ). (A) Absorption image of the atomic cloud 
after probing the cavity for 0.5 ms and subsequent ballistic 
expansion over 7 ms. (B) Excited atomic population N e with 
momenta (±hk, izhk), deduced from the regions enclosed by 
blue circles in A, as a function of the probe-pump detuning 
5 (circles). The solid line shows a fit, based on the theo- 
retical model (see text). The background in N e originates 
from a residual thermal component extending into the de- 
tection region. (C) Recorded (circles) and fitted (solid line) 
mean intracavity photon number n p h during probing. The 
dashed line indicates the intracavity photon number of the 
probe pulse, obtained in the absence of atoms. (D) Temporal 
evolution of N B (t) inferred theoretically from the fit shown in 
C. (E) Mean intracavity photon number (n p h), averaged over 
the probe interval, as a function of 8 (circles). The fit (solid 
line) takes into account the calibrated photon number of the 
probe pulse (dashed line). From the recorded photon signal 
we deduce the photon number originating from Bragg scat- 
tering off the created excitations (dashed-dotted line). The 
shading in B and E indicates the fluctuations caused by vari- 
ations of the relative phase ip between different experimental 
runs. The detuning of the pump beam from the empty cavity 
resonance was A c = — 2n x 18 MHz and the total atom num- 
ber N = 1.6(2) x 10 J . The damping constant 7 was set to 
0.6kHz (see text). 
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where & creates particles in mode |e). The first term in 
Eq. [2] corresponds to the kinetic energy of atoms in mode 
|e) , with the recoil energy E T — \£r and the atomic mass 

771. 

The elementary excitations of the Hamiltonian H are 
of collective nature and their eigenenergy E s softens for 
increasing interaction strength V (Fig. IB). The excita- 
tion energy E s vanishes at a critical interaction strength 
V CT , where interaction and kinetic energy are balanced, 
i.e. iV|T4r| = 2E T . This marks the transition point be- 
tween a normal and a supersolid phase. In the supersolid 



phase, a macroscopic occupation of the mode |e) gives rise 
to a checkerboard pattern in the atomic density distribu- 
tion, which is accompanied by the build-up of a coherent 
cavity field amplitude an. The energy of collective ex- 
citations rises again with increasing interaction strength 
and approaches the single-particle excitation energy of 
the induced optical checkerboard potential. 

To measure the excitation spectrum of the system at 
momenta (±hk,±hk), we perform a variant of Bragg 
spectroscopy (24j [26]. After preparing the system at a 
given interaction strength, the cavity field is excited with 
a weak probe pulse along the x-axis (Fig. 1A). Interfer- 
ence between the cavity probe field and the transverse 



pump field results in an amplitude-modulated lattice 
potential r]^n pr cos(dt + tp) cos(kx) cos(fcz), with probe- 
pump detuning S = u> pl — u> p , relative phase ip and mean 
intracavity photon number n pl of the probe field. In 
momentum space, this corresponds to the perturbation 
H pr = rjy / n pT Ncos(St + <p)($ + c). After probing, all 
laser fields are switched off. This projects the created col- 
lective excitations onto the free-space momentum states 
\±hk,±hk), which are detected via absorption imaging 
subsequent to ballistic expansion (Fig. 2A). A typical 
resonance curve of the excited momentum state popu- 
lation iV c = (c^c) as a function of S is shown in Fig. 2B. 
Clear resonances are revealed both for positive and neg- 
ative probe-pump detuning, corresponding to stimulated 
scattering of probe photons into the pump field and vice 
versa. The corresponding resonance frequency E s /K is 
obtained from a fit (solid line) based on the model de- 
scription Ea. [2] [25]. 

The measured excitation spectrum as a function of 
pump power P is displayed in Fig. 3. When increasing 
the interaction strength in the normal phase towards the 
critical point, the excitation energy E s exhibits a dis- 
tinct softening. In contrast, for positive V, the excita- 
tion gap is observed to increase with interaction strength 
in accordance with the absence of a phase transition. 
The increasing gap reflects the tendency of the cavity- 
mediated interactions to suppress density fluctuations 
with A-periodicity. 

We monitor the cavity output light on a single-photon 
counting module. This gives us real-time access to the 
formation of collective excitations during probing. Be- 
cause of matter-wave interference between the ground- 
state component of the condensate and the created ex- 
citations, a spatial density modulation evolves in time. 
Transverse pump light is Bragg scattered by this den- 
sity modulation into the cavity mode. The output from 
the cavity therefore follows the oscillatory evolution of 
this density modulation (Fig. 2C). The total intracav- 
ity photon number, averaged over the probe interval, 
again exhibits a double resonance (Fig. 2E) , whose char- 
acteristic shape has its origin in interference between the 
probe field (dashed) and the Bragg scattered pump field 
(dashed-dotted) . This provides a second method to de- 
termine the resonance frequency E a /h, as shown in Fig. 3. 

We model the dynamics of the system during prob- 
ing by evolving the Hamiltonian H + H pl in time, tak- 
ing into account the recorded shape of the probe pulse, 
see Fig. 2C. Quantitative agreement with the data is ob- 
tained by including a phenomenological damping rate 7 
in the time evolution of c |25j . The corresponding time 
evolution of the momentum population N c (t) is shown 
in Fig. 2D. Depending on the relative phase ip, which is 
not controlled in the experiment, the phase of the excited 
density oscillation varies between different experimental 
runs. This leads to an intrinsic fluctuation of the quan- 
tities N e and (n p h) , as indicated by the shaded areas in 
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FIG. 3. Excitation spectrum. Measured resonance frequen- 
cies E s /h, obtained from atomic {N e ) and photonic ({n p h}) 
signals, are shown in blue and red, respectively, for positive 
(open circles) and negative (filled circles) interaction strength 
V. Gray shading shows the theoretical prediction including 
experimental uncertainties 25 . The experimental parameters 
are N = 1.7(2) x 10 5 , A c = -2tt x (19.8, 23.2)MHz for V < 
(normal and supersolid phase) and A c = +2-7T x 15.1 MHz for 
V > 0. For V > 0, the pump power is scaled in terms of the 
critical pump power observed for A c = — 2-w x 19.7 MHz. 



Fig. 2B and E. 

In the supersolid phase, the perturbation creates exci- 
tations on top of the macroscopic steady-state population 
in momentum mode |e) [32] . Because of matter- wave 
interference, the detected population N c is strongly af- 
fected by fluctuations of the relative phase <p. Yet, the 
variance of N c around the steady-state value, deduced 
from several experimental runs, shows clear resonances. 
Alternatively, we use the oscillation amplitude of the in- 
tracavity photon number around the steady-state value 
|an| 2 to extract the excitation energy E s 25 . As shown 
in Fig. 3, the excitation energies deduced from these two 
signals rise again with pump power P and provide a con- 
sistent picture. 

The observed excitation spectrum (Fig. 3) is in qual- 
itative agreement with a soft mode energy E a — 
2E r y/l + V/\V CI \ which follows from Eq. V2\ Quantita- 
tive agreement is obtained by further taking into account 
the lattice potential of the transverse pump beam as well 
as contact interactions between colliding atoms [25] . Ac- 
cordingly, the single-particle mode |e) in Eq.^ns replaced 
by a Bogoliubov mode |1) and the interaction energy 
gets renormalized by the matrix element of V(r, r') be- 
tween the condensate mode and the excited state |1). 
In the normal phase, the state |1) lies in the lowest en- 
ergy band of the transverse lattice potential with quasi- 
momentum (±hk,±hk) and bare energy E\. For P = 
the excitation energy E s reaches the Bogoliubov energy 
E\ = h x 8.7 kHz, which is shifted by the mean-field en- 
ergy with respect to the kinetic energy 2E,- = h X 7.5 kHz. 
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FIG. 4. Response to density perturbations. Shown is the 
density response in the normal and supersolid phase, normal- 
ized to the non-interacting case (P = 0). The shaded areas 
show the theoretical prediction including experimental uncer- 
tainties and fluctuations caused by the uncontrolled relative 
phase if. The atomic damping constant 7 = 2tt x 0.5(1) kHz 
was obtained from the fit to the data in the normal phase. 
The inset shows the response of the excited momentum state 
population iV e for negative (filled circles) and positive (open 
circles) interaction strength V. The experimental parameters 
were the same as in Fig. 3. 



phase, in accordance with Feynman's relation. In the 
supersolid phase, the density response decreases again as 
a function of P, as expected for the increasing excitation 
energy (Fig. 3). A similar behavior is observed for the 
response in the momentum state population iV e (Fig. 4, 
inset), which exhibits a larger sensitivity to variations 
of the relative phase <p [25] . For positive interaction V 
(Fig. 4, inset), the response of the system is reduced with 
respect to the non-interacting case, V — 0. 

We have observed a roton-type mode softening causing 
a superfluid-to-supersolid transition in a model system 
for long-range interactions. Increasingly complex spatial 
structures of long-range atom-atom interactions [28H30] 
can be tailored by extending the experimental setup to 
multiple cavity modes. 
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In the supersolid phase, the Bogoliubov mode |1), which 
is dominantly affected by the long-range interactions, lies 
in a higher energy band at zero quasi-momentum of the 
emerging checkerboard lattice potential. As a result of 
competition between the increasing band energy and the 
renormalized interaction energy, the excitation energy E s 
rises again. 

Another signature for a softened dispersion relation at 
finite momentum is provided by the susceptibility to ex- 
ternal density perturbations. In terms of the static struc- 
ture factor S (k), which measures the density-density re- 
sponse, this is expressed by Feynman's relation S(k) = 
(Hk 2 /2m)/uj(k), where w{k) denotes the dispersion of a 
homogeneous system |27j . To determine the density re- 
sponse in our system, we quantify the amount of pump 
light which is Bragg scattered into the cavity. This al- 
lows us to measure the probe-induced density modulation 
(Sp e ), where Sp c denotes the density fluctuation operator 
at momenta (±hk,±hk). We define the corresponding 
quadratic density response function 1Z P as the spectral 
weight of (Sp c ) 2 , normalized to the pump power P and 
the probe pulse area (n pr ). The response function IZn 
of the atomic population N c in the excited momentum 
state |e) is defined in a similar way |25j . 

Approaching the critical point from below, a strongly 
increasing density response is observed (Fig. 4). This 
indicates the presence of enhanced density fluctuations 
with A-periodic correlations in the unperturbed system. 
The data is in good agreement with our model (shading) , 
which predicts 1Z P to scale as (Ei/E s ) 2 in the normal 
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SUPPLEMENTARY MATERIALS 

We give details of the experimental setup and pro- 
vide a theoretical description of the dispersively cou- 
pled condensate-cavity system. An effective Hamiltonian 
for the cavity-mediated long-range interaction is derived. 
Based on a mean-field description, we numerically cal- 
culate the steady state of the system including contact 
interactions and the transverse lattice potential. We de- 
duce in a Bogoliubov approach the energy of collective 
excitations, which softens towards the critical point. The 
time evolution of the system during probing is derived 
and a diverging response to density perturbations at the 
phase transition is found. 

Experimental details 

The atoms are prepared in the hyperfine state 
(F,rriF) — (1,-1) with respect to a quantization axis 
pointing along the cavity axis, where F is the total angu- 
lar momentum and mp the magnetic quantum number. 
The transverse pump laser with wavelength A = 784.5 nm 
is linearly polarized along the y-axis (see Fig. 1 in the 
main text) and off-resonantly drives, via the atoms, two 
degenerate cavity TEMoo modes with circular polariza- 
tions e + and e_. The ratio of the corresponding two- 
photon Rabi frequencies is given by ??+/??- = 3.26/1.25, 
where all allowed dipole transitions in the D\ and D 2 
lines of 87 Rb have been taken into account. The max- 
imum dispersive shift of the two cavity modes induced 
by a single maximally coupled atom is Uq = 2tt x 87 Hz 
and Uq = 2ir x 33 Hz, respectively. The transverse pump 
laser induces an optical lattice potential along the z-axis 
with periodicity of A/2. Its depth V p is calibrated using 
Raman-Nath diffraction [3T] and takes for our experi- 
mental parameters a value of 3£" r at the critical point, 
with recoil energy E T = %^-, wavevector k — 2n/\, and 
atomic mass m. In the theoretical analysis, the Gaussian 
envelopes of the pump and cavity fields along the trans- 
verse directions [35] are effectively taken into account by 
weighted averages of V p , rj± and U over the spatial ex- 
tent of the atomic cloud. The Thomas-Fermi radii of the 
condensate in the external harmonic trapping potential 
are given by (R X ,R V ,R Z ) = (3.5, 18.3,3.7) um, assuming 
an atom number of N = 1.65 x 10 5 . 

The length (176 /xm) of the cavity is actively stabilized 
using a laser with a wavelength of 830 nm, which is refer- 
enced onto the transverse pump laser. The depth of the 
resulting intracavity lattice potential was measured to 
be 0.04(2)i? r , and is neglected in the theoretical analysis. 
Transverse pump light and cavity probe light propagate 
through independent optical fibers, resulting in a varia- 
tion of their relative phase ip between different experi- 
mental runs. The cavity output light is monitored on a 
single-photon counting module with an overall detection 



efficiency of intracavity photons of 4(1)%. For the data 
taken in the normal phase, the critical pump power P CT 
was deduced from the intracavity photon number moni- 
tored during independent sweeps across the phase tran- 
sition. 



Theoretical description of the coupled 
condensate-cavity system 

The coupled condensate- cavity system is described in a 
many-body formalism following Ref. '33, 341- By adiabat- 
ically eliminating the fast cavity field dynamics, we derive 
an effective Hamiltonian, which describes the long-range 
atom-atom interaction. 

After adiabatically eliminating the electronically ex- 
cited states, the transversally driven condensate-cavity 
system is described by the many-body Hamiltonian H = 
H c + H a + Ha,- C , where 






HA c a*a 
d 3 r¥{r) 



f* f (r)*(r) 



2m 
*(r 



+ V P cos 2 (kz) 



(SI) 



H„ 



d 3 rW (r)[hrj cos(kx) cos(kz)(a + a/) 
ft[/ cos 2 (fcr)a t a]\['(r) , 



with bosonic atomic field operator $(r) and photon oper- 
ators a and a*. To keep the notation simple we describe 
only one of the two circularly polarized TEMoo cavity 
modes. Final results will be given for the case of two 
cavity modes. 



In equation (SI I, H c describes the dynamics of a single 
TEMoo cavity mode with spatial mode function cos(fcx), 
whose frequency uj c is detuned by A c = u; p — lu c from 
the pump laser frequency w p . The term H a captures the 
atomic evolution in the transverse optical lattice poten- 
tial with depth V p , including contact interactions with 
strength q — 47rfi a , where a denotes the s-wave scatter- 
ing length. The interaction between the atoms and the 
pump and cavity light fields is governed by H a - C . The 
first term describes light scattering between pump and 
cavity field at a rate which is determined by the max- 
imum two-photon Rabi frequency rj. The second term 
accounts for the dispersive shift of the cavity resonance 
frequency with light-shift U$ of a single maximally cou- 
pled atom. 

As the cavity field reaches a steady-state on a time 
scale fast compared to atomic motion, its equation of 
motion can be formally solved, yielding 



(Ac 



776 
U B) 



(S2) 



with the cavity field decay rate k = 2tt x 1.25 MHz. 
Due to Bragg scattering of pump light, the intracav- 
ity field amplitude is proportional to the order param- 
eter 9 = J d 3 r ^(r) cos(fca;) cos(kz)^>(r) which mea- 
sures the atomic density modulation on the checkerboard 
pattern cos(fcx) cos(kz). The overall dispersive shift of 
the empty cavity resonance caused by the presence of 
the atoms is proportional to the bunching parameter 
B = J d 3 r &t( r ) cos 2 (kx)^{r). 

After eliminating the steady-state cavity field of 
Eq. (p2J) from Hamiltonian Eq. (SI), an effective Hamil- 



tonian description is obtained (see main text, Eq. 1): 
fleff = H a + V J d 3 rd 3 r'*t ( r )i£rt ( r ') C os(kx) cos(fcz) 

x cos(fcr') cos(/cz')^(r)^(r'). 

(S3) 
The strength V of this cavity-meditated atom-atom in- 

2 A 2 

teraction is given by V = h-S- — %■ « h^—, where the 

J A 2 +k 2 A c ' 

detuning of the pump laser from the dispersivcly shifted 
cavity resonance A c = A c — U0B0 was taken to be large 
compared to the cavity half-linewidth k. Here, Bq = (B) 
denotes the bunching parameter in the steady state. In 
a quantized picture, NV/h corresponds to the rate at 
which cavity photons are exchanged between atoms, as 
shown exemplarily in the zoom of Fig. 1A in the main 
text. 



The effective Hamiltonian Eq. (S3) describes a closed 



system and neglects dynamical and quantum backaction 
effects originating from cavity decay and cavity input 
noise. This is justified on short timescales as long as 
|A C | > « » £ r /7i [35]. 



Mean-field description in the steady state 



Based on a mean-field description [36], we derive a 
numerical solution of the steady state of the system. 

A mean-field description of the system is obtained by 
formally replacing the operators \fr and a in Eq. (SI) 
with the atomic mean-field 



^Nipo and the coherent cav- 
ity amplitude ao, respectively. Their steady-state values 
are determined by the non-local Gross-Pitaevskii equa- 
tion 



H ip (x,z) 



2m 



(d 2 x + d 2 z ) + V p (z)+hU (x)\a f 



+ thri(x,z)(ao + a Q )+g2T>\ipo\ )^Po(x,z) 



with ao 



V&o 



(S4) 

and the chemical potential /j,q. Here, 
we introduced the steady state order parameter 0o = 
iV(^ | cos(fca;) cos(A;,z)|^o) an d bunching parameter B = 
N(ip \ cos 2 (kx)\ipo) , as well as the notations V p (z) — 
V p cos 2 (kz), r](x,z) = ?]cos(kx)cos(kz) and Uq(x) = 



Uq cos (kx) . We reduced the description in Eq. ( S4 ) to 



the pump and cavity directions, assuming a homogeneous 
condensate density along the weakly confined y-axis. The 
contact interaction strength is accordingly replaced by 
<?2D = A 2 n<? with average 3D condensate density n and 
the normalization condition L dx L dz \ipo\ 2 = 1 [3T] . 



For negative A c , Eq. (|S4j) exhibits a dynamical insta- 
bility above a critical transverse pump power P cr , driving 
the system from a normal phase into a supersolid phase 
with A-pcriodic density modulation. In the normal phase, 
P < P C r, the condensate density is flat along the cavity 
axis, resulting in a vanishing order parameter, 0o = 0, 
and cavity field amplitude, ao = 0. The mean-field solu- 
tion ijjq is given by the lowest energy Bloch state in the 
shallow optical lattice potential of the transverse pump 
field. In the supersolid phase, P > P cr , the atomic cloud 
exhibits a A-periodic density modulation both along the 
transverse and the cavity direction. Correspondingly, the 
order parameter takes a finite value, Oo 7^ 0, and light 
scattering off the diagonal Bragg planes results in a co- 
herent cavity field amplitude, a =/= 0. The mean-field 
solution ifio is given by the minimal energy state in the 
two-dimensional lattice potential originating from inter- 
ference between the transverse and cavity field. In the su- 
persolid phase, the emergent checkerboard density mod- 
ulation locks to one of two possible sublattices, which 
are spatially shifted by A/2 and have opposite signs of 
Go and cto. As both solutions share the same excitation 
spectrum, we assume Go > in the following. 

We numerically find the ground state ipo of the sys- 



tem by propagating Eq. (S4| in imaginary time, using 
a computational cell of size A 2 with periodic boundary 
conditions. We include the recorded steady state mean 
intracavity photon number |ao| 2 and the experimentally 
calibrated lattice depth V p in the calculations. 



Deriving the collective excitation spectrum 

We describe collective excitations on top of the mean- 
field solution. To this end, we calculate the Bogoliubov 
excitations in the steady-state lattice potential and iden- 
tify a single excitation mode which is dominantly affected 
by the long-range interactions. We diagonalize the re- 
sulting truncated Hamiltonian and deduce the energy of 
collective excitations. 

Following Ref. |36| . we expand the atomic and cavity 
field operators around their mean-field solution (ipo,ao) 
according to 

* = (VNi/jo + 5^)e" ittlo/h 

a — ao + 5a. (S5) 

Here, the linearized cavity field fluctuation operator 



5a = 



V 



A r 



se 



GoU 



-SB 



(S6) 



is given in terms of the operators SO = y/N J d 3 r (S 1 ^^ + 
8$) cos(kx) cos(kz)t[) Q and SB = y/Nfd 3 r(5& + 
8^f) cos 2 (kx)tp , describing fluctuations of the order and 
bunching parameter around their steady-state values Go 
and Bq, taking -0o to be real- valued. 

After eliminating the cavity field fluctuation Sa, as- 
suming A 2 3> K 2 , and keeping only linear terms in the 



equation of motion for 
Hamiltonian 



<5 'J, we arrive at a quadratic 



H,. 



H + V 



SO 



QqUq 
A, 



SB 



(S7) 



describing the dynamics of atomic fluctuations around 
the mean-field solution. The first term Hq captures the 
dynamics in the static steady-state lattice potential in- 
cluding contact interactions, and is given by 

h2 -{dl+dl) + V v {z) 



Pn = / d 3 r S& 



2m 



hr)(x, z)(ao 



HU {x)\a \ 



<5* 



>V>o( 



92B% [SV + VV) 



h2 



2g 2U \^ Q \ 2 S^5i>. 
(S8) 



The second term in Eq. (S7| describes how atomic den 



sity fluctuations are affected by the cavity-mediated long- 
range interaction. 

In the normal phase, H cxc is directly obtained from 



Eq. (S3) by quadratic expansion in (5^ around the steady 



state. In the supersolid phase, the dispersive cavity shift, 
which depends on the atomic density distribution, gives 
rise to an additional long-range interaction term with 
A/2-periodicity along the cavity axis [31], whose strength 
is proportional to the order parameter 0o- However, as 
the order parameter vanishes at the critical point, solely 



<5*(r) = J2 ( u i( r ) £ i + v *A r )4) > 



the term SO 2 in Eq. ( S7 ) induces the normal to supersolid 
phase transition. 

To find the energy of collective excitations described 
by Pexc, we first calculate the elementary excitations of 
H using the Bogoliubov ansatz [38l [39] 



(S9) 



with band index j, Bogoliubov modes ttj(r) and Vj(r), 
and corresponding mode operators Cj, fulfilling bosonic 
commutation relations. In this basis the Hamiltonian 
P exc reads 

Pcxc = ]T (E^c, + NV Xj (Cj + cj) 2 ) , (S10) 



with Bogoliubov energies Ej and interaction matrix ele- 
ments 



Xj = (V'o 



cos(kx) cos(fcz) + 



e p ( 



cos 2 (/ex) 



Uj + vj 



(Sll) 
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J2n=± VnXn, with V± and x± deduced accordingly from 
the cavity parameters {jj± , U ) . The calculated energy 
spectrum shown in Fig. 3 takes into account the exper- 
imentally calibrated depth of the transverse lattice po- 
tential and the measured steady-state intracavity photon 
number |ao| 2 in the presence of atoms. Systematic un- 
certainties of these quantities are estimated to be 10% 
and 25%, respectively (see shaded regions in Fig. 3). 



FIG. 5. (A) Bare energies (blue) of the excited state with 
maximum matrix element \ in the normal and supersolid 
phase. Dashed and solid blue lines correspond to the case of 
vanishing and present contact atom-atom interactions. The 
red dashed-dotted line shows the soft mode energy in the pres- 
ence of long-range and contact atom-atom interactions. (B) 
Maximum matrix element \ in the presence (solid) and ab- 
sence (dashed) of contact atom-atom interactions. 



From a numerical calculation of the matrix elements Xj 
we find in our parameter regime only a single Bogoli- 
ubov mode exhibiting a relevant matrix element Xj f° r 
energies up to 15E T . All other matrix elements Xj are 
suppressed by more than two orders of magnitude. In 
the normal phase the maximally coupled excited mode 
lies in the lowest energy band of the transverse lattice 
potential with quasi-momentum (q x ,q z ) — (±hk,±hk). 
In the supersolid phase this state is folded into the third 
band at the center of the Brillouin zone as a result of the 
emerging A-periodicity. 

In terms of the excited mode with dominant matrix 
element Xj > which we denote in the following by the index 
j = 1, the effective Hamiltonian reduces to 



H r 



E 1 c\c 1 + NV X (ci + c\) 



ha 



(S12) 



with x — Xi- I n t ne absence of the transverse lattice 
potential and contact atom-atom interactions, H cxc re- 
duces in the normal phase to the Hamiltonian given in 
the main text, Eq. 2. 

We obtain the soft mode energy spectrum by diagonal- 
izing Hamiltonian H exc in terms of a second Bogoliubov 



Probing the collective excitation spectrum 

By integrating the equation of motion of the atomic 
field, we predict the build-up of atomic excitations and of 
the cavity light field during probing. 

To probe the excitation spectrum, we weakly drive 
the cavity field with amplitude r/ pr (t) and frequency 
Lu pl . This is described by the driving Hamiltonian 
-hn pr {t)(ae^ St+ ^ + a)e-^ St+ ^) with probe-pump de- 
tuning 5 = w pr — Up and relative phase ip between probe 
and pump beam. The resulting coherent intracavity 
probe field is given by 



a pT (t) = - 



Vpr(t)t 



-i(St+ip) 



A r 



(S14) 



Interference of the probe field with the transverse pump 
field and the steady-state cavity field cto results in a mod- 
ulated lattice potential. The corresponding perturbation 
of the atomic field is given by 



H n 



Hit) 



66 



UqQq 
A r + ih 



SB 



cos(6t + ip), (S15) 



with perturbation amplitude £(£) = 2rjy/n pr (t) and mean 



intracavity photon number n pT (t) 



As |A C | > k 



in the experiment, we set the phase shift, originating from 
coupling into the cavity, to n. In terms of the mode 
operator ci, the perturbation reads 



H pl = h^(t)y/Nx(ci + c[)cos(St + ip) 



(S16) 



transformation b 



(J,Ci 



st 



Up to a constant term, the Hamiltonian H = H P 



To quantify the response on this perturbation, we evolve 



this yields H exc — E s Wb with soft mode energy 



-ffp r in time according to 



E s 



£W1 



ANVx 



(S13) 



the Hciscnbcrg equation 

ihci =E x Cx + 2NVx(c! + c\) 

+ H{t)^/Nxcos{8t + <p)- ifvycx 



(S17) 



For V < 0, E s softens towards the phase transition 
and vanishes at V cr = —E"/(4:Nxci) with (£?"', x cr ) = 
(E% , x) I P=P x ■ I n t ne supersolid phase, the excitation en- 
ergy rises again due to competition between the increas- 
ing Bogoliubov energy E\ and the decreasing interaction 
energy which is proportional to the matrix element x ( see 
Fig. %. 

Results for two circularly polarized cavity modes a± 



Motivated by our experimental observations, we phe- 
nomenologically introduced a damping term with damp- 
ing constant 7 into the time evolution of c\. This ac- 
counts for possible damping mechanisms like s-wave scat- 
tering with other momentum modes, trap loss or finite- 



size dephasing. The general solution of Eq. (S17) reads 



are obtained by replacing Vx hi (S12) and (S13l by 



ci( 



E 



(t) = 2 V ^N X n pT ,o( -±Im(y(t))+iRe(y(t)) . (S18) 



E. 
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Here, 



y(t) = e (^-T)* / dt' e'^-^' cos(St' + p)U(t') . 



(S19) 

with u! s = E s /h and n pT (t) — n pr ,oTl(t) 2 . Here, the en- 
velope function Ti(t) of the probe pulse with duration r 
and the maximum probe photon number n prj o were intro- 
duced. The population of the excited momentum mode 
and the mean intracavity photon number are directly ob- 



tained from Eq. (|S18|) 



(&\<k) 



4r] 2 n prfi Nx 



El 



lm(y(t)) 2 +Re(y(t)f 



n p h(t) 



a 



4»7 \f'i 



A 



hrfi N X (E X 



IK \E; 



Im(y(t)) 

2 



(S20) 



The second term in n p h(£) describes the pump field which 
was Bragg scattered off the excited density modulation 
into the cavity mode. From N\, the population N c in 
the momentum state |e) is obtained according to iV c = 
(Ni , where ( denotes the absolute square of the Fourier 
amplitude of the excited Bogoliubov mode u\ + v\ at 
momenta (±hk,±hk). 

The curves shown in Fig. 2C and D of the main 
text are obtained from Eqs. (pi), where the phase ip 
was adjusted to fit the data in Fig. 2C. The resonance 
curves, displayed in Fig. 2B and E of the main text, 
correspond to phase-averaged values (A ir c (r)) ¥ , e [ .2 7 r] an d 
(J din p h(i)/r) ve [o i 27r] where the resonance frequency 
w s and the amplitude rj 2 Nx were adjusted independently 
to fit the data. Due to the detection background an offset 
was added to N c (t). The standard deviation of fluctua- 
tions associated with the uncontrolled relative phase ip is 
displayed by the shadings in Fig. 2B and E of the main 
text. 



Response in atomic density and momentum state 
population 

Bragg spectroscopy, as described by the density pertur- 
bation H pr , induces a response in the atomic density and 
in the population of the excited state. Normalized to the 
integrated amplitude of the perturbation, this density re- 
sponse is a measure for the susceptibility of the system 
on an external density perturbation and provides at zero 
temperature a direct link to the static structure factor 
\38l \40j - We use the recorded modulation (So) of the 
intracavity field to quantify the density response ((5/5) of 
the system. 

Consider the density operator p(r) = ^^(r)^'(r) and 
its linear expansion p — po + Sp = N\ip \ 2 -\- VlVV'o (<5^ + 



5^) around the equilibrium value po = (p) . In terms of 
the density fluctuation operator in Fourier space, defined 

as SfJk — T^S^^-ikrAA 

and bunching parameter 60 and SB read 



Jd 3 re lkr <5/5(r), the fluctuations of the order 



se = Yl <W 4 

ke(±fc,±fc) 



and 



5B= J2 <W 4 - 

ke(±2fc,0) 



Fluctuations of the cavity field 5a, see Eq. (S6), are thus 
given by 



da = 



4(A C 



ke(±fc,±fc) 



e i/ 

A r + in 



ke(±2fc,0) 



and provide a measure of the induced atomic density 
modulation. We define the corresponding (quadratic) 
density response function 



K p = 



1 



Pfin pI (t)dt 



dS 



\(Sa)\ 2 dt 

« |7,/(4(A c +iK))P 



ve[0,27r] 



where again an average over the relative phase p is per- 
formed. 

Experimentally, we extract |((5a)| 2 from the resonance 
fits based on Eq. (J3J) (see Fig. 2E in the main text), 
with experimentally calibrated |«o| 2 and n pr . In order to 
probe the linear response of the system, the probe power 
was accordingly lowered when approaching the critical 
point. The lowest perturbation applied in the normal 
phase corresponds to 30(8) intracavity probe photons in 
total. 

Similarly, the response function IZjy associated with 
the detected number of atoms with momenta (±hk, ±hk) 
is defined as 



n N = 



l 



PJodtn pT (t) 



d5(iV e (r)) ye [o,2 W ]- 



which is proportional to the area below the fitted reso- 
nances of N e (r), see Fig. 2B in the main text. 

The dominant scaling factor of 1Z P when approaching 
the critical point is found from Eq. pi) to be (Ei/E s ) 2 . 
As 1Z P measures the quadratic density response, this is 
in agreement with Feynman's relation, given in the main 
text. Since N c (t) is sensitive to both quadratures of c x 
(see Eq. ([3])), the response function 1Z N scales differently 
and exhibits larger variances in the relative phase p>. 
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